Conservation of spin current 
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The conventional definition of spin-current, namely spin density multiplied by the group velocity, 
is not a conserved quantity due to possible spin rotations caused by spin-orbit (SO) interaction. 
However, in a model with spin-spin interactions, rotation of a spin causes a dynamic response of 
surrounding spins that opposes the rotation. Such a many-body effect restores the spin-current 
conservation. Here we prove that the non-conservation problem of spin-current can be resolved if a 
self-consistent spin-spin interaction is included in the analysis. We further derive a spin-conductance 
formula which partitions spin-current into different leads of a multi-lead conductor. 
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Recently, considerable interest has been paid to the 
quantum physics of spin-curren1ji. It is believed that 
a controlled spin-current generation, detection, and us- 
age can provide interesting applications to spintronics. 
Spin-current generation has been classified as "extrin- 
sic" or "intrinsic" . An extrinsic spin-current is gen- 
erated by external physical factors and driving forces 
of the spintronics device, such as optical spin injection 
achieved experimentally and the various spin pumps 
studied theoretically^. An intrinsic spin-current is gen- 
erated by physical factors existed inside the spintronic 
device, notable is that generated by various spin-orbit 
interactions^*SiL&. In particular, it has been theoret- 
ically predicted that non-magnetic systems with spin- 
orbit interaction and under an external electric field, can 
generate a spin-current flowing perpendicular to the elec- 
tric field£. Such an spin-current is termed "dissipation- 
less" because the electron motion is perpendicular to the 
electric field. There are so far extensive theoretical work 
on spin-current physic al 10 , and some experimental works 
have also appearec&Ai which may provide support to 
some of the theories. 

Despite the increasing literature on spin-current 
physics, it is recognized that the definition of spin-current 
itself is still somewhat controversial. If one mimics the 
definition of charge-current, then a spin-current 7 S can 
be defined as the time derivative of spin density N s , 
I s ^< dN s /dt >, here the bracket is the quantum aver- 
age. In its simplest form at steady state, such a definition 
gives a spin-current I s = |(Jf — 7;), where Jf, /j. are the 
charge current with spin-up and spin-down, respectively. 
Clearly, this is a very intuitive definition of spin-current 
and is adopted by most of the work in literature. Since 
both spin and velocity are vectors, the spin-current is a 
tensor. In systems where there is a spin-orbit interaction, 
the spin density is not conserved: spins can rotate from 
their initial orientation due to the interaction. Therefore, 
spin-current 7 S becomes a non-conservative quantity. A 
quantity which is not conserved is difficult to study ex- 
perimentally and indeed, it is unclear what is even mea- 
sured if an experimental detection method can be found 



to measure spin-current. Without knowing what is mea- 
sured, the definition of spin-current becomes non-unique 
and there have been several definitions in literature. 

Although it is unclear a priori if a measurable spin- 
current must be conserved, the property of conservation 
would be nice to have, at least theoretically. Such an 
issue has been discussed in a recent paper of Sun and 
Xiei^ who defined an additional rotational contribution 
to spin-current on top of the above conventional defini- 
tion, and the total spin-current is then conserved. How- 
ever, there is no microscopic theory concerning this prob- 
lem. In the present work, we further investigate this 
problem by using the conventional definition of spin- 
current 7 S ^< dN s /dt >, but we include a spin-spin in- 
teraction into the Hamiltonian for computing the quan- 
tum average < • • • >. Our basic idea is the following. 
The existing problem of spin-current non-conservation 
was due to SO interaction which rotates the spin away 
from its original direction. In a spin system with spin- 
spin interactions, the rotation of one spin causes a dy- 
namic response of surrounding spins through the inter- 
action, and this response opposes the original spin ro- 
tation. This many-body effect introduces a new term 
in the spin-current that balances the spin-current non- 
conservation. In the end, one obtains a final spin-current 
in the interacting model (which has a different value than 
that of non-interacting model) that is conserved. In the 
following, we demonstrate this idea by proving that 7 S 
obtained in a model with dipole-dipole interaction is, in- 
deed, conserved. Importantly, the idea is based on a 
many-body phenomenon, as long as there is spin-spin 
interaction, our conclusion should be general. Further- 
more, since there are always some dipole interactions be- 
tween spins, including their contribution is also rather 
natural. Finally, the conserved spin-current gives a lin- 
ear spin-conductance, we derive an expression for this 
spin-conductance for multi-probe device systems. 

We consider a coherent multi-probe mesoscopic device 
which consists a device scattering region that connects 
to a number of leads extending to far away. In a lead la- 
beled by a, the spin-current I sa ^< dN s /dt >. Following 
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Refii^, we calculate the time derivative using Heisenberg 
equation of motion and in steady state, the spin current 
I s can be written in terms of Green's functions: 

I** = -/ ^Trh(G^<-E<G° 

+ G<£«-5^G<)] (1) 

where S r,a are the retarded and advanced self-energies 
due to the presence of leads; G <,r '° are the lesser, re- 
tarded and advanced Green's functions of the device scat- 
tering region, respectively. Note that the trace in the last 
equation can be reduced as: 

Tr [cr3(G r S< - £<G a )] = Tr[G a FG' r S<] 

where 

F= [GT^s-^GT 1 = Kff ]-SV 3 -a 3 57 . 

Here Hq is the Hamiltonian of the device scattering re- 
gion without the leads. Hence Eq.QJ can be written as 

Is* = - J ^Tr[G<(£> 3 -^) 

- G<CE a a 3 -a 3 S r ) + [a 3 ,H }G<] (2) 

where G< = iGT Q / Q G a with T@ the linewidth func- 
tion. 

The total spin-current is 

= = (3) 

a 

where 

dS . [dE_ ^ 

~dt ^ 1 J ~4n ■ ^ 

with 03 = — i[er3, H ]. Eq.© is just the continuity equa- 
tion for spin current in the absence of external magnetic 
field and spin relaxation^. The physics of this equation 
is rather clear: in the presence of any spin-orbit interac- 
tion the electron spin precesses due to an internal torque 
on the spin, hence the total spin current flowing into the 
scattering region is non-zero and it equals to the rate of 
spin precession —dS/dt. This is as if there were some 
"spin accumulation" in the scattering region. In other 
words, the total spin-current I s is itself not conserved 
due to spin precessing. 

This situation is reminiscent to the continuity equation 
of charge current I e in the presence of a time dependent 
field 15 , namely 

'-£'- = -f ( 5 ) 

a 

Here the total charge current I e is equal to the charge 
accumulation in the scattering region. The problem of 



charge current conservation under AC fields was dis- 
cussed by Buttikcr 1 - who pointed out that the total par- 
ticle current is not conserved under AC conditions, but 
the total particle current plus total displacement current 
is a conserved quantity. The total displacement current 
is precisely dQ/dt in the above equation. Since dis- 
placement current results from induction which is related 
to electron-electron interaction, Biittiker formulated a 
current conserving theorjii^ by including the electron- 
electron interaction at the mean field level, which natu- 
rally deduces the displacement current in each lead. The 
problem of displacement current partition in multi-probe 
conductor under noncquilibrium conditions has been re- 
ported in Refill within the nonequilibrium Green's func- 
tion formalism. 

The similarity between Eqs.10 and suggests that 
the problem of spin-current conservation can also be 
looked at from an interacting spin point of view. As dis- 
cussed above, a spin-spin interaction causes a dynamic 
response inthe system which counter-react on any rota- 
tion of a spin. In the following, we demonstrate this 
many-body physics by investigating the consequence of a 
self-consistent spin-spin interaction in a device Hamilto- 
nian which contains the Rashba spin-orbit (SO) interac- 
tion. Indeed, we prove that the self-consistent spin-spin 
interaction produces a term that exactly cancels the right 
hand side of Eq.@ so that Isa = 0. 

Our model Hamiltonian is (ft = 1): 

H = H + V ss 

where 

H = £(e„ + qU n )dld n + ^(C44 + h.c.) (6) 



H " = 2 £ ( V nLA d m4dj + h.C.) 
nmij 

where H ss is the spin-spin interaction and n,m,i, and j 
include spin indices. Here V so and V ss , respectively, are 
matrix elements of the Rashba SO interaction and the 
spin-spin interaction. In real space, they are given by 

V so (x)=a R a-(zxk) (7) 

andi£ 

V ss (x lX ') = {gii B ?o ■ V a ■ V . 

\x — x I 

To deal with transport problems, we make a mean field 
analysis on the spin-spin interaction so that H ss becomes 

H ss = V nm did m 

n m 

with 

Vn m = \ £i » < d\d, > +h.C.) 
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In real space, V nm becomes, 

V(x) = gi-i B <y ■ V0m = -gvB<? ■ H 

with 



>M 



V F - x'\ 



(8) 



(9) 



The local magnetic moment M is determined by spin 
density which, in the language of nonequilibrium Green's 
functions (NEGF), is given by 

/dE 
— Tr s [aG<J 

where Tr s is the trace over spin space and G% x is the 
diagonal element of G < . Note that H = —V4>m is the 
self-consistent effective magnetic field due to the spin- 
spin interaction. The Green's function also depends on 
this effective field through V(x), 



G r = 



1 



E-Hq- g[i B a ■ Vcj) 



M 



(10) 



From Eq. (JHJ , 4>m satisfies the following Poisson-like equa- 
tion 



V <t> M = -4%pM 



(11) 



where pj^ = — V • M is the effective density of magnetic 
charge^. These equations form a self-consistent problem. 

To solve Eg. we consider a very large volume of 
space surrounding the device scattering region such that 
the total magnetic charge inside that volume to be zero. 
Mathematically, this consideration means: 

JdxV-M(x)=0 or J ^Tr(V -crG^) = 0(12) 

where Tr = Tr s Tr Q includes the trace over both spin 
space and orbit space. Ea.i[T2[l can be achieved if the 
volume in which we solve Ea. fll|> is so large that the 
effective magnetic field H on the surface of that volume 
is zero. 

Using the model Hamiltonian Eq. JBJ and applying the 
Heisenberg equation of motion, we can evaluate 
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-i[<Ta,H ] = 2ia R y ■ a 



This result together with Eq.l|12[) proves 

f -o- 

Therefore the "spin accumulation" of Eq.QJ is actually 
zero. In other words, if the spin-spin interaction is in- 
cluded in the Green's function, there will be no spin accu- 
mulation and the spin current is conserved. We empha- 
size that in the self-consistent formalism, Eas. (jlll()lll|l 
form the basic set of equations for the spin current con- 
serving theory. 



Now we derive a "spin conductance" Q that corre- 
sponds to the spin-current. When the external bias is 
small, we expand the spin-spin interaction V(x) in terms 
of bias v a , V(x) — ^2 a u a v a where we have introduced 
the notion of spin dependent characteristic potentiaU^i^ 
u a which satisfy the gauge invariant condition ^2 a u a — 
1, i.e., the spin current depends only on the difference in 
external bias. Expanding Eq. (J2J in terms of small bias 
Vp, we find at zero temperature, 



where 



Isa = ^ Ga/3 Vp 



Qafi = Tr[g Q/3 - C^2g ai )xxU(j{x)] 



(13) 



(14) 



where the matrix g a p is given by 



[S a cr 3 — a 3 'E r )5 a p 
[a 3 ,H o ]S a0 ]G r T p G a + h.c. (15) 



Note that the first term of Eq. I|14|) comes from expansion 
of Fermi distribution function and the second term in- 
volving the characteristic potential is due to the spin-spin 
interaction term in the expansion. Linearizing Ea. (|lll) wc 
have 



- v 2 4> a {x) = kTi- s [(> y g V j)o 



c u a (x) - (y^g T a)xx](16) 



where k = (2-KgiiB)l olr and <f>M = ^}2 a 4>aV a - Note 
that the expansion over external bias, the spin conduc- 
tance g a p in both Ea. (|15(l and Ea. (|16(l do not depend 
on the self-consistent interaction. From Eq.© we have 
u a = g/J-BP • V<^>a- From Poisson like equation Ea. i|ltj|) 
we obtain the spin dependent characteristic potential and 
the spin conductance can be calculated from Eq. i|14|) . 

We emphasize that this conductance guarantees that 
the linear spin-current of Ea. (|13fl is conserved. With- 
out the spin-spin interaction, the conductance would be 
given by only the first term on the right hand side of 
Eq. (|17|) , and the resulting spin-current would not be con- 
served. In fact, the sum of the second term of Eq. l|17|) 
over space, i.e., the quantity Tr[^ 7 g 7 J, is exactly equal 
to the "spin accumulation" —dS a /dt in the scattering re- 
gion from lead a due to Rashba interaction in the small 
bias limit. To see this, we find from Ea. i|15|) 



dS a 
dt 



-Tr[^2g ia ] = ^Tr[h,i/o]GT Q G"] 
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so that ^ a dS a /dt v a = dS/dt. This means that the 
spin-spin interaction puts this contribution into the spin 
conductance itself automatically, so that the right hand 
side of Eq.© vanishes. Eq. l|17|l therefore partitions the 
non-conserving part of the spin-current (i.e. the right 
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hand side of Eq.JJJ) when there is no spin-spin interac- 
tion) into each leads, such that the spin-current becomes 
conserved. 

Since the solution of the Poisson like equation requires 
numerical calculation, analytically we can avoid this by 
using a quasi- neutrality approximation, i.e., assuming 
that the effective density of magnetic charge pm{%) = 
so that the local magnetic moment M is independent of 
position. Then, we find the spin dependent characteris- 
tic potential by setting the right hand side of Eq. (|16[) to 
zero and obtain: u a (x) = (£) g ja )xx/{J2m 9r n )xx- The 
conductance is then found to be: 



Gap = Tr[g a/3 



7 9ja) 
(Si ( 7 5i)7) 



(17) 



This expression is very similar to that of charge current 
partition in ac situationsi^ii, and it partitions the to- 
tal spin-current into each lead a so that the total spin- 
current flowing into the device is conserved. 

The above microscopic theory result, Eq. l|17|l . is valid 
for Rashba SO interaction. For a general SO interac- 
tion, a similar expression to Eq. (|17Jl can be derived using 
a phenomenological argument^. To do that, we require 
two conditions: (i) the total spin current is conserved; (ii) 
the value of spin current depends only on the difference of 
external bias. The latter condition means that spin cur- 
rent remains unchanged if external bias at each lead is 
shifted by the same amount. Now, the unconserved spin 
current I° a is given by Eq. i|13|) . The "spin accumulation" 
If ee dS/dt is given by £ Q I c sa = T,p(E a G c a0 )vp = -if, 
where we have used G%n = Trg a p for the non-conserved 
spin conductance. Note that the total "spin accumula- 
tion" is due to the contribution lf a from each lead a, i.e. 
I~ = J2 a Ifa- Since only the total "spin accumulation" is 
known, we need to find lf a by partition the spin current. 
For this purpose, let I sa = I% a + A a lf, or equivalently 



where A a is an unknown to be determined. Condition 
(i) gives J^a Gap = 0, hence we obtain J2 a ^a — 1. Con- 
dition (ii) gives gauge invariance Yp Gap = 0, hence we 
obtain Tr[J^ g a p — A a J2jp 9ip\ = from which we find 
A n = Y J1 Qo, 1 /Y Jlv t 3iv Therefore Eq.JEJ gives: 



GaP = G a f 



(Y G c 



(19) 



Ea. (|19() has the same form as Ea. H17|) which is specific 
to Rashba SO. We therefore propose that Eq. l|17fl can 
serve as a phenomenological theory which conserves spin- 
current regardless of the detailed spin-orbit interactions. 
Namely, if we use Ea. (|17f) to compute spin-conductance, 
the resulting spin-current from Eq. l|13l) will always be 
conserved regardless of which SO interaction is present. 

In summary, we have proven that the conventional 
spin-current I s ^< n s v s > for Rashba interaction be- 
comes a conserved quantity if spin-spin dipole interac- 
tion is included. Such a dipole interaction introduces 
a self-consistent field which correlates spins spatially. 
For general SO interactions, a phenomenological theory 
for spin current partition is proposed which conserves 
spin current, and the resulting spin-conductance has the 
same form as that derived from the microscopic theory 
of Rashba interaction. 
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